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Abstract 

Using nuclear magnetic resonance (NMR) techniques with three-qubit sam- 
ple, we have experimentally implemented the highly structured algorithm for 
the 1-SAT problem proposed by Hogg. A simplified temporal averaging proce- 
dure was employed to the three-qubit spin pseudo-pure state. The algorithm 
was completed with only a single evaluation of structure of the problem and 
the solutions were found with probability 100%, which outperform both un- 
structured quantum and the best classical search algorithm. 
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I. INTRODUCTION 



Quantum Computers exploiting the principles of quantum parallelism and entanglement 
of quantum states outperform classical counterparts greatly. Quantum algorithms can solve 
computational problems much faster than the classical ones |l|-f§], and can solve even some 
intractable problems to classical methods, such as the most famous problem of factoring 
large integers [[|. Hitherto, the selected systems to realize quantum computation have been 
proposed such as trapped ions quantum dots 0, and cavity quantum-electrodynamics 
MM- However, the most effective one to realize quantum algorithms to date is the liquid- 



state NMR spins Several quantum algorithms have been implemented by liquid- 

state NMR ensemble including the Deutsch-Jozsa algorithm with two, three [|T^ -[l8[| and five 
|19| qubits, the Grover's algorithm with two p0|j2l[] and three qubits fl22 |, an algorithmic 



benchmark ||23|| , an order-finding algorithm with five qubits [ [24|j . 

Search algorithm is one of the most important algorithms. The difference between quan- 
tum and classical search algorithms is that, quantum ones can operate simultaneously on a 
superposition of all possible search states so as to be superior to classical counterparts greatly. 
A quantum search algorithm that uses the underlying structure of the search problem to 
establish correlations between problem properties and the solutions results in a substantial 



improvement over a search one that ignores them, such as Hogg's algorithm |25[] vs. Grover's 
one Q. However, when only two qubits are considered, there is no difference between Hogg's 
algorithm and Grover's one because both of the algorithms only require one search step. For 
n > 2 qubits, the superiority of Hogg's algorithm will be displayed, which still require a sin- 
gle step, while Grover's one requires v^2™ steps. For example, when n = 3 qubits, Grover's 
algorithm needs at least 2 steps with probability near 95%. Hogg's algorithm for 2-qubit has 
been implemented in Ref. |[26|| . In this paper, we experimentally implemented the Hogg's 



algorithm for a 1-SAT with the 3-qubit pseudo-pure state prepared by three-step temporal 
averaging. 
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II. THE HIGHLY STRUCTURED ALGORITHM 



The highly structured algorithm is associated with the satisfiability problem (SAT) ||27|| , 
one of the most difficult class of non-deterministic polynomial problem (NP) f2~7| . A SAT 



is a combinatorial search problem J27J, consisting of a logic formula in n boolean variables 
(true or false) Vi } ...,V n and the requirement to find an assignment, specifying a value for 
each boolean variable, that makes the formula true. The logic formula can be expressed as 
a conjunction of m clauses and each clause is a disjunction of some variables, any of which 
may be negated. When all the clauses have k variables, the problem is called fc-SAT. In 
general, the computational cost of solving a SAT grows exponentially with n in the worst 
case. For k ^ 3, fc-SAT is NP-complete. However, for 1-SAT and maximally constrained 
k-SAT, Hogg proposed a quantum algorithm to find a solution in one step by using problem 
structure. As each false clause for a given assignment is counted as a conflict, solutions are 
assignments with no conflict. Therefore, the correlations between problem properties and 
the assignments are set up so that Hogg's algorithm can solve problems with only a single 
testing of all the assignments, corresponding to a single classical search step. 

For a 1-SAT problem with n variables and m clauses, in fact, the algorithm takes the 
following four stages [p5[] . 

1. Initialize the amplitude equally among the states, giving an equal superposition of 
bases \s), = 2~% J2 S \s) , with bit strings s being all 2 n assignments of n variables. 

2. Adjust the phases of \s) based on the conflicts c in the assignments s, ranging from 
to m, i.e., apply a transformation R on \ipi), where R is a diagonal matrix with elements 

{y/2 cos[(2c — for even m; 

1 4J (1) 

i c , for odd m. 

3. Mix the amplitudes from different assignments with the mixing matrix U depending 
only on the Hamming distance d between r and s is described as 

. 2 2~ cos[(n — m + 1 — 2d)|], for even m; 

U rs = U d (r,s) = { 1 (2) 

2 -t e i7T(n-m)/4(_^d ) for odd m. 



Here, d(r,s) = \r\ + \s\ — 2\r A s\ with \r\ (|s|) being the number of 1-bits in r(s) and 
\r A s\ , the number of 1-bits both assignments have in common. The operator U can be 
defined in terms of two simpler operations: U = WTW, where W = 2~5(— l)l rAs is the 
Walsh-Hadamard transform, , and T is a diagonal matrix with elements 

{\/2cos[(m — 2h — for even m; 

(3) 

depending only on the number of 1-bits in each assignments h = \r \ , ranging from to m. 
4. Measure the final superposition \ipf) = UR\ipi). 

For our 3-spin system (i. e. n = 3), the number of clauses , m, in the 1-SAT formula can 
be 3, 2, or 1. The quantum circuit for the highly structure search in a 3-qubit system is shown 
in Fig. 1. According to four stages of the Hogg's algorithm above, the calculation starts 
with a Walsh-Hadamard transform W applying to a pseudo-pure state 1 000) , to prepare the 
equal superposition state \ipi), where W = Ha <8> H B ® if<7,with H = X 2 Y (pulses applied 
from right to left) being a single-spin Hadamard gate. Then we derived the NMR pulse 
sequences for R and U in Eqs .(1) and (2). We take the 1-SAT with m — 3, i.e., the logic 
formula being V\ A V% A V3, as an example. In this case, U could be represented by WTW. 
From Eqs.(l) and (3), the expressions for R and T can be written into 

diag [R Vi av 2 av 3 } = H, -1, -1, h -1, h h 1], (4) 
diag[T m=3 ] = [1, i, i, -1, i, -1, -1, -i]. (5) 

which correspond to the pulse sequences Rv^^/Wz ^ ZiZ 2 Z 3 , r m=3 ^ ZiZ 2 Z% (up to an 
irrelevant overall phase factor). In order to make the best use of the available coherence time 
and to diminish errors due to the increased number of RF pulses, the pulse sequences of W, 
R and F were optimized to eliminate unnecessary operations with the help of NMR principle 
|28| . Two instances for m — 1 and 3 were chosen for demonstrating Hogg's algorithm in our 
experiments for their simple pulse 

sequences. All the reduced pulse sequences for all possible formulas with m = 1 and 3 
as well as theoretical solutions are listed in Table. 1. 
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III. THE NMR EXPERIMENTAL RESULTS 



The Hogg's algorithm was implemented by liquid-state high-resolution NMR spec- 
troscopy with carbon-13 labeled alanine NH^ — C a il(C^i2 2 ) — C Oj dissolved in _D 2 0. 
We chose C , C a and C 13 as the three-spin system in experiments, representing spin 1, spin 
2 and spin 3, respectively. Spectra were recorded on a Bruker ARX500 spectrometer with 
a probe tuned at 125.77MHz for l3 C. The chemical shifts of three different carbon spins 
were calibrated about —4320, 0, and 15793Hz, and the spin-spin coupling constants Ji 2 , J23, 
and J13, 34.94, 53.81, and 1.21Hz, respectively. The relaxation times were measured to be 
11(1) = 20.3 sec, 21(2) = 2.8 sec, and 21(3) = 1.5 sec and T 2 (l) = 1.3 sec, T 2 (2) = 0.41 sec, 
and T 2 (l) = 0.81 sec fL8| . Protons were decoupled during the whole experiments, using a 
standard heteronuclear decoupling technique. Spin-selective excitations were executed using 
low-power, long-duration pulses with a Gaussian shape. The length of these pulses was tai- 
lored to achieve sufficient selectivity in the frequency domain without disturbing the nearest 
nucleus, depending on the difference of the chemical shifts between nuclei. We took the 
length of the selective pulses to be 0.7ms and the excitation power to be 29.6dB for selective 
I pulses for all 13 C nuclei. 



A. The preparation of pseudo-pure states 

In liquid-state NMR ensemble quantum computers, Chuang et al. |TT| proposed pseudo- 
pure states instead of the genuine pure states as the initial state. Several methods have 
been proposed to prepare the pseudo-pure states including spatial averaging [JlOfl , temporal 



averaging (T^|, logical labeling [|Il1 , ^0|, p9l . Recently, an alternative simplified temporal aver- 
aging was proposed by Lieven et al. ||24|| . For a n-qubits homonuclear system, the deviation 
density matrix of the system in thermal equilibrium is a sum of n terms, 



n 



Peq = ^2hz- (6) 



k=l 



The deviation density matrix of the desired pseudo-pure state |00...0) is 



n n n 

P e ff = ^ kz + Ikzllz + Ikzhzlmz + , (7) 

fc=l «>fc=l m>l>k=l 

a sum of 2 n — 1 terms. In order to obtain the desired pseudo-pure state, [(2 n — l)/n] different 
experiments are at least needed by utilizing controlled-NOT {CN^ flips spin j if and only 
if i is |1) ) and NOT (iVj flips the sign of spin i). 

For a homonuclear 3-spin system, (i.e. n = 3), three separate experiments are at least 
required to achieve the pseudo-pure state p oo- We use 3 separate experiments, giving a 
total of 9 product operator terms. The 2 extra of 9 terms are canceled out pairwise by 
using NOT (iVj) operations. The procedure can be chosen in a variety of ways. For different 
samples, one can choose an optimal experimental scheme. As to our sample of alanine with 
1/2 Ji 3 = 0.41 sec, which is comparable to the smallest T 2 , gate CN\ 3 or CN31 containing 
the evolution under the scalar coupling J13 should be excluded in the preparation circuit. 
We choose the following 3 states preparation sequences acting on the equilibrium as 

E, CN 32 CN 21 N 3 , CN 21 CN 12 CN 32 (8) 

where E is unit matrix denoting no operation. The corresponding product operators of the 
3 states are 

Pi = hz + hz + hz 

P2 = ^hzhzhz + Zhzhz - hz ( 9 ) 
P3 = Zhzhz + 2/12/22 + hz 

The sum of the 3 states is 

Psum = Pi + P2 + p3 = hz + hz + hz + ^hzhz + ^hzhz + ^hzhz + ^hzhzhz (10) 

It can be seen from Eq. (7) and Eq. (10) that p sum is the same as p e ff for n = 3. The 
experimental pulse sequences was shown in Fig. 1. Of course, the same result to exclude 
this 1/2 Ji 3 evolution can also be gained by swapping gates proposed by Collins et al. [ITS 



but it is more complex to implement experimentally. Note that a magnetic field gradient 
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pulses were accomplished at the end of the experiments in order to eliminate the residual 
transverse magnetic vectors and thereby reduce the errors. 

The pseudo-pure state pooo was achieved in the experiment, by summing the three ex- 
perimental results, where the qubits from high to low corresponds to spin 1, 2 and 3, respec- 
tively. The resultant spectra of the three spectra then recorded (shown in Fig. 2) through 
reading-out pulses confirm that a pseudo-pure state pooo has been prepared. The normalized 
diagonal elements of the deviation matrix p oo for the pseudo-pure state 1 000) by quantum 
state tomography []30| was given as 



diag[p 000 ] = [1.000, 0.0314, -0.0291, -0.0032, 0.0520, 0.0114, -0.0535, -0.0277] (11) 

All diagonal elements except one with the value of 1 in Eq. (11) should be zero theoretically. 
So the maximal relative error of the experimental values of the diagonal elements was shown 
to be < 6%, with the small off-diagonal elements. 

Similarly, the procedure can be applied to a homonuclear 4-spin system (i. e. n = 4.) 
To obtain the pseudo-pure state poooo, we can perform 5 experiments as follows 

CN 12 CN U CN 31 , CN 21 CN A2 CN U} CN 12 CN A2 

(12) 

CN 12 CN U N 3 , CN 23 CN M N 31 

Among a total of 20 terms produced, the 4 extra terms can be eliminated pairwise through 
Ni operations and one extra term removed through the magnetic field gradient technique. 
For a homonuclear 5-spin system, the procedure has clearly been demonstrated in Ref. |f24| . 



B. The NMR results of the highly structure algorithm 

Applying the pulse sequences in Table. I. to the pseudo-pure state p 000 , we got the 
results of Hogg's algorithm. For m = 3, because only the diagonal elements are non-zeros, 
we reconstructed the normalized diagonal elements like the pooo preparation above. The 
experimental deviation matrixes were obtained respectively 
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diag[p VlAV2AVs \ = [-0.0190, 0.0297, -0.0582, 0.0631, -0.0072, 0.0416, -0.0800, 1.0000], 
diag[pv lAV2 AV 3 ] = [0.0087, -0.0074, 0.0959, -0.0845, 0.0105, -0.0056, 1.0000, -0.0393], 
diag[p VlA v 2AV3 \ = [-0.0412, 0.0716, -0.0149, 0.0187, -0.0580, 1.0000, -0.0047, 0.0289], 
diag[p VlAV2A y 3 ) = [0.0037, 0.0340, -0.0186, 1.0000, -0.0211, 0.0092, -0.0670, 0.0881], 

(13) 

diag[p VlAV2AV3 ] = [0.0512, 0.0077, -0.157, -0.0269, 1.0000, -0.0127, 0.0029, -0.0082], 
diag[pv lA v 2 AV 3 ] = [0.0173, 0.0171, 1.0000, -0.0348, -0.0091, -0.0092, 0.0606, -0.0093], 
diag[p VlA v 2A v 3 ] = [-0.0304, 1.0000, -0.0197, 0.0200, -0.0228, 0.0594, -0.0199, 0.0199], 
diag[pv lA v 2 AV 3 ] = [1-0000, 0.0314, -0.0291, -0.0032, 0.0520, 0.0114, -0.0535, -0.0277]. 

In contrast with the theoretical expectation, the maximal relative errors of the experi- 
mental values of the diagonal elements was shown to be < 9%, with the small off-diagonal 
elements. 

For m = 1, we reconstructed the experimental final deviation matrices by quantum state 
tomography . In Fig. 4 were shown the theoretical and measured results for these cases 
with the maximal relative errors 16% ~ 25% of the density matrix elements. It can be seen 
clearly from Eq. (13) and Fig. 4 that the experimental results are in good agreement with 
the theory. Errors are primarily due to inhomogeneity of RF fields and static magnetic fields, 
magnetization decay during the measurement and imperfect calibration of the rotations. 

IV. DISCUSSION 

In summary, we have experimentally demonstrated a NMR realization of quantum algo- 
rithm for a highly structured searching problem on a three-qubit quantum computer. The 
simplified temporal averaging for preparing the pseudo-pure states takes less separate ex- 



periments than the original temporal averaging proposed by Knill et al. [ I3"|| , e. g., 3 times 
instead of 7 times for a 3-spin system. Moreover, its Signal-to-Noise ratio is better than 
that of spatial averaging [[K],|r^]. Hogg's algorithm for 3 qubits was completed on the pre- 



pared pseudo-pure state pooo, unlike some algorithms only on the thermal state, and the 
final results were also read out by a weak measurement on the ensemble, which accomplish 
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the whole process of quantum computation. 

Hogg's algorithm based on the structure of a problem, with a particular choice of the 
phases determined by the number of conflicts in assignments, finds the solutions to 1-SAT 
problem in a single step with probability 100% even for n — > oo . By contrast, the un- 
structured search methods require 0(2 n ) steps classically, and 0(2 n l 2 ) steps on quantum 
computers @]. For example, in the 3-qubit Grover's algorithm NMR experiment imple- 
mented by Vandersypen et al. |22| with a heteronuclear system, per Grover iteration applied 



19 pulses, 2 evolutions of 1/2 J and 3 of 1/4 J and performed 2 iterations at least to find 
a single solution \x$). Comparing our pulse sequences with it, Hogg's algorithm adopts far 
less RF pulses and requires far less coherence time so as to be realized more easily. In 
addition, probability 100% guarantees that Hogg's algorithm is a complete method, i. e. 
failure to find a solution definitely indicates the problem is not soluble. Unlike the previ- 
ously proposed quantum algorithms, such as the Grover's algorithm ||, that find solutions 
with probability less than one, cannot guarantee no solutions existing. Hogg's algorithm 
is more efficient to solve any 1-SAT problem. Furthermore, this algorithm also applied to 
maximally constrained soluble fc-SAT problems [25| for any k with an analogue procedure. 



To experimentally implement this algorithm for larger systems, the main difficulties are to 
address and control the qubits well and to maintain coherence during the computational 
process. 
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Captions of the figures 



Fig. 1 Quantum circuit for 3-qubits that implements a highly structure search algorithm. 
Horizontal lines represent qubits, time goes from left to right. Using three Hadamard gates, 
an pseudo-pure states \ipo) = 1 000) is transformed into a uniform superposition state 
which is then converted to the answer state \ipf > after the action of gates R and U. For 
the definition of R and U, see text. 

Fig. 2 NMR pulse sequences to implement the pseudo-pure state p oo- Narrow and wide 
boxes correspond to f and ir pulses (refocusing pulses) respectively. X and Y denote the 
pulses along the x- and y-axis, X and Y, opposite to the x- and y-axis. (a) to prepare the 
state of p2 in Eq. (4). (b) to prepare the state of ps in Eq. (4). p\ is the thermal equilibrated 
state, no operation). The pulse sequences are designed for alanine with small J 13 . 

Fig. 3 The resultant 13 C spectra of three experimentally measure spectra of the pre- 
pared pseudo-pure state p oo- The reading-out pulses (a) [| ] (b) [|]^ , (c) [|]^ were applied 
respectively. The abscissa indicates the frequency, and the ordinate denotes the intensity of 
the spectra (in arbitrary unit). 

Fig. 4 Experimental and theoretical deviation density matrices (in arbitrary units) for 
the Hogg's algorithm of all logic formulas when m — 1. (a) — (f) represent the recovered 
matrices for the logic formula Vi, Vi, V 2: V 2 , V 3: V 3 , respectively, the left and right column 
denoting the real and imaginary components, (al) — (fl) are the corresponding theoretical 
values. Relative errors are shown as percentages. 
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Captions of the tables 



Table. I All the logic formulas for 1-SAT of a 3-spin system for m = 1 or 3 , the 
corresponding theoretical solutions and the reduced pulse sequences for the operator U RW, 
where the logic variables with subscript i stand for spin i. X and Y denote the | pulses 
along the x- and y-axis, X and Y, opposite to the x- and y-axis. The subscripts represent 
the qubits. 
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Fig. 1 xhpeng et al. 
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Fig .2 xhpeng et al. 
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Fig. 3 xhpeng et al. 
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